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WELL-POSEDNESS OF SPACE AND TIME DEPENDENT 
TRANSPORT EQUATIONS ON A NETWORK 


ARNE ROGGENSACK 


Abstract. This article is concerned with the study of weak solutions of a linear 
transport equation on a bounded domain with coupled boundary data for general 
non smooth space and time dependent velocity fields. The existence of solutions, 
its uniqueness and the continuous dependence of the solution on the initial and 
boundary data as well on the velocity is proven. The results are based on the 
renormalization property. At the end, the theory is shown to be applicable to the 
continuity equation on a network. 


1. Introduction 


In this article, we study linear transport equations on bounded domains with cou¬ 
pled boundary values and prove their well-posedness. As a main application we have 
in mind transport phenomena on a network which occur in various fluid dynamical ap¬ 
plications. In that case, on each edge of a graph a linear transport equation is imposed. 
The equations are coupled at the nodes by linear boundary conditions mapping from 
the outflow part of the boundary to the inflow part. From an abstract point of view, 
one can see this as a multivariate transport equation on a one-dimensional domain 
with coupled boundary data 


Pt + (Up), + C P = f 


in (0, T) x 

in Q 
on Tt* 


p(0, ■) = Po 

{v\J)-p={v\J)-U{p\r T ) 


We will prove existence, boundedness, uniqueness, and stability of solutions p G 
U([0,T],L p ((0, l)) n ) for non-smooth matrix-valued velocities U = diag(w) with u G 
^((OjT), 1U 1 ’ 1 ((0, l)) n ) and u x G L 1 ((0, T), L°°((0, l)) ra ) without any restrictions on 
the sign of the velocity. In the Lipschitz continuous setting, every change of the sign 
of the velocity implies a change of the direction of a characteristic. This can lead to 
additional difficulties since it has two effects: 

First, even if the data is smooth, there does not have to exist a classical smooth 
solution as the domain is bounded. For example, by the method of characteristics the 
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solution of 


on (0,1) x (0,1) is given by 


pt + m -1 ) P ) X = o 


p(t,x) 


Po (X - t(t - 1)) 

Pin -x'j 


if t < \ or x > (t — i) 2 
otherwise, 


which has possibly a discontinuity if no further conditions are imposed. The only way 
to get smooth solutions is to require additional conditions on the initial and boundary 
data po and p; n , which depend on the characteristics and thus on the velocity. In many 
applications however, the velocity is not known a priori. 

The second issue concerns the coupled boundary data. In this case, we can compute 
the characteristics, but due to the coupling conditions, we are not able to use them 
directly to find a solution if the velocity has an infinite number of changes of the sign 
in a finite time. On a network, this can lead to a characteristic of the form of an 
infinite tree. Each time the characteristic intersects an inner node it divides itself into 
multiple parts. 

For these reasons we will use a different and very general approach looking for weak 
solutions. This approach is based on the concept of the renormalization property, 
which was introduced in 1989 by Lions and DiPerna J] for tangential velocity fields. 
This concept was extended to bounded domains in with inflow boundary conditions 
and velocity fields with a kind of Sobolev regularity by Boyer [2] in 2005 and Boyer 
and Fabrie [3]. Recently, these results were generalized to velocity fields with BV 
regularity by Crippa et al. and to stochastic differential equations by Neves and 
Olivera m- For constant velocity fields, the transport equation on a network was 
studied by Sikolya in her dissertation [IS] in 2004. In 2008, this was expanded to 
infinite networks by Dorn mm- Especially, the long time behaviour was analysed by 
a semigroup approach. 

In this article, we generalize the results of Boyer and Fabrie to the case of coupled 
boundary data. In particular, we will study a vector valued transport equation with 
affine linear coupled boundary conditions and a time and space dependent velocity 
field. 

Our results for the transport equation can e.g. be applied to study the existence 
of solutions of a low Mach number model on a network. Details to this topic can be 
found in a forthcoming paper or in the PhD thesis m- However, transport equations 
on networks occur also in various other contexts as traffic flow H2HU, gas flow in 
pipe lines PQ or river flow mm 

The article is divided in six sections. In Section [2J we introduce the notations and 
requirements. Section[3]recalls the renormalization property and it provides uniqueness 
and boundedness results. In section 31 we prove the continuous dependence of the 
solution on the data before we show the existence of a solution in Section [5j In the 
last section, we illustrate how these general results can be applied to the continuity 
equation on a network. 
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2. Assumptions and requirements 

In this section, we introduce the precise setting. For simplicity, we will only consider 
the one-dimensional domain 52 = (0,1) and denote the boundary T = <952 = {0,1} and 
r T = [0, T] x r. Since we have in mind to consider one-dimensional edges of a graph 
this is no restriction for the application. The outer normal vector on T is called v and is 
given by v = — (—l) w for oj G T. We use the integral notation f r although the boundary 
consists only of two single points, i.e. d ui is the counting measure. An advantage of 
this notation is the easy generalization to higher dimensions. All considerations in this 
article could also be done for the case of bounded multidimensional domains 52 with 
Lipschitz boundaries. 

The aim of this article is to find and characterize a solution p(t, x) G R" of the 
problem 

p t + (Up)* + C p = f in (0, T) x 52 

p(0, •) = po in 52 (1) 

(uU)-p=(vU)-n(p\r T ) onT T 

with initial conditions po G A 00 (52)" with po > 0 almost everywhere and with an affine 
linear boundary operator T-L. We assume 

uGL 1 ((0,T),W 1 > 1 (52)"), 

c € T 1 ((0, T) x 52)", 

{u x + c)~ GL 1 ((0,T),L oo (52)"), (2) 

{u x ) + € T 1 ((0,T),L oo (52)") 

and 

/ € T 1 ((0,T),L°°(52)") 

with / > 0 and denote by U = diag(u) and C = diag(c) diagonal matrices with the 
diagonal entries u 3 and c 3 , respectively. On Ty we define component-by-component 
the measure 


dp u — (uu)dwdt (3) 

and introduce its positive part d //+ = (uu) + dwd2, its negative part dp~ = (isu)~dujdt 
and its absolute value |d/i„| = d/x+ + dp~ . For each j, these measures divide the 
boundary Ty into two parts, the inflow part with uuj < 0 and the outflow part 
with vuj > 0. For p G [l,oo) the space L p (Tt, dp,±) = L p {Tt, R", dpj) is provided 
with the norm 

\\g\\p,w,± = J ^(|. 9 |) T W(u'u) ± dwd 2 

with the function /3(x) = (x\ • • • ai(() T and with a positive definite diagonal weight 

matrix W G R" x ". Here, the absolute value has to be understood component-by¬ 
component. We do also consider the space L°°(Tt , d p,±) = L°°(Tt , R", dp±) equipped 
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with the maximum norm of the componentwise L°°-norm, i.e. 

/ 


|oo,± = max 
j 


ess sup I I 

. («,t)er T 

\ u 3 '(t, w) ± 5^0 


The boundary operator % has to assign a boundary value to the inflow part of the 
boundary, i.e. T~L is a mapping 

%: L°° (Tr, dp. J) —>■ L 00 (Tt , d/4“). 

We assume the mapping to be affine linear, i.e. it is given as H(p) = pi n + Q{p) where 
it is pi n £ L°°(r T ,dp-) with p ln (t,x) > 0 for dp“-almost all ( t,u ). Here, Q is a linear 
operator fulfilling the following conditions: 

- The operator Q : L°°(Tt , dp+) —>• L°°(rT, d p~) is weakly-* continuous, i.e. for 
all p n £ i°°(TT,dp+) with p n p in L°°{TT,dp +) it holds 


Q{Pn) ^ Q{p) 


in L°° (T t , dp u ). _ 

- There is a positive definite diagonal weight matrix W £ R nxn such that the 
L 1 -operator norm of Q is less or equal one, i.e. for all p £ L°°(IY, dp+) it holds 


\\G(P)\\ 


1 ,w,— 




- ^Tlie operator Q is causal, i.e. for all p £ L°°(Tt , dp+) and almost all t £ [0, T] 
it holds 

X[o,t]G(X{o ,t]P) = X[o ,t]G{p). (4) 

- The operator Q is dp+-almost everywhere positive, i.e. it holds Q{p) > 0 dp“- 
almost everywhere for all p £ L°° (Tt, dp+) with p > 0 dp+-almost everywhere. 

- There exists a constant vector p max £ K™ such that component-by-component 
the inequalities 


Pin exp 


a(s)ds 


F+ + G{F+) + G (Pmax) < Pmax 


(5) 


for dp u -almost all (t,w) £ Tt and 


PO — /An ax 


for almost all x £ fl are valid. Here, it is a(s) = ||(u a; (s, •) + c(s, •)) ||ioo(Q)n 
and 

/||(/i(s,-)) + lk~\ 

ds. 

As we will see later, this construction of T~L is motivated by the transport equation on 
a network. 


F. 


.(t) = J exp J a(r)drj 
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To conclude this section, we observe three properties of the operator Q. First, for 
all p £ L°°(rT,dpu) an d all t £ [0 ,T] the inequality 


r \g(p)\ T w(isu)~du J dt= [ [ X[0A \g( P )\ T w(vu)-dujdt 

/o Jr Jo Jr 

= [ f |^(x[o,t]P)| T W(uu)"dwdt 

Jo Jr 

= ll^?(X[0,i]P)l|l,iu,- 
< ||X[0,i]Pl|l,tu,+ 
r T r 

IX[o,t]P| T W(uu) + dwdt 


ia Jr 

r t 


|p| T W(uu) + du;df 


( 6 ) 


o Jr 


is true. Second, the equation Q holds in fact for all essential bounded functions 
g £ L°°([0,T]): Let p € L°°(rT, d/z+), p / 0 and £ > 0 be given. Since the step 
functions are dense in L°°([0, T]) we can choose a step function g n = Ylk=i a kX[t k ,t k+1 ] 
with 


llff - gn IIl°°([o,t]) < 


e 

2||p||i,uj,+ 


Then, because of g{f n p ) = /raf?(p) it follows 


II g{gp) - gg(p)h,w- < || g{{g - g n )p)\\i,w- + ll(ff - g n )G(p) lli,™,- 

< II (s - 9n)p\\l,w,+ + || g - 5n||L“([0,T])ll<7(p)l|l,™ 

< \\g - 9n\\L°°([0,T])\\p\\l,w,+ + II g - ffra||L“([0,T])||p||l,uj,+ 

< £ 


and thus 

g{gp) = gg{p ) (?) 

holds d/x“-almost everywhere. 

The third observation concerns the positivity of the operator. Using the positivity 
of g and the triangle inequality, we estimate for p £ L^Tt, d/x+) 

ig(p)) + = \{g{p) + \g(p)\) 

= \ (g{ P + ) - g(p~) + 1 g(p + ) g ( P ~)|) 

<\{g{p + )-g{p-) + \g(p + )\ + \g{p-)\) 

= g(p + ) 

almost everywhere in Ty. In the same way, we can prove that 


(G(p))~ < G(p~) 


(9) 
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holds. In the following, we often extend functions in L p (Ty, d/xj) by zero to con¬ 
struct functions in L p (Tt, |d/z„[) or L p (ry ) n ■ Especially, we extend, if necessary, the 
operators Q and TL to operators mapping from L°°(Tt , |d// u |) —>• L°°{Tt , |d/z u |). 


3. Uniqueness and Boundedness 

In this section, we will prove the uniqueness and boundedness of solutions of the 
initial-boundary-value problem CD- To this end, we first introduce the concept of 
weak solutions before we discuss regularity and the renormalization property of the 
solutions. Equipped with this, we are able to prove uniqueness and boundedness. For 
the latter, we will give explicit lower and upper bounds. 

We start with the usual definition for a weak solution. 


Definition 1 . We call a function 


p£L°°{[0,T) x!l)" 

a weak solution of the transport equation p t + (U p) x + C p = / if it holds 

0=1 [ P T (<p t + U(/? x - Cip) + f T ipdxdt 

Jo J n 

for all ip £ C 0,1 ([0, T ] x Cl, KL”) with <^(0, •) = ip(T, •) = 0 and tp = 0 on Ty. 


A priori, it is not clear how to define boundary conditions for this class of weak 
solutions since p is not a Sobolev function and thus, the usual trace theory is not 
applicable. 

In [3], Boyer has proven some properties of such weak solutions for scalar equations, 
especially the existence of a trace and the renormalization property of the solution (see 
the following theorem). 

Theorem 2 (Trace theorem, Boyer 3j). Let Cl C R d be a d-dimensional bounded 
Lipschitz domain and let be u £ L 1 ((0,T),W 1 ’ 1 (Cl) d ) and c, / £ L 1 ((0,T) x Cl) with 
(c T u x )~ £ L 1 ((0, T), L°°(f2)) and ( u x ) + € L 1 ((0, T), L°°{Cl)). Then, for each weak 
solution p £ L°°((0,T) x Cl) of the scalar transport equation 

p t + di v(up) + cp = f 


the following properties hold: 

Time continuity: The function p lies in C([0, T\,L p (Cl)) for all p £ [1, oo). 
Existence and uniqueness of a trace function: There exists a unique essen¬ 
tially bounded function 'yp £ L°°(Yt, |d/i K |), called trace, such that for any 
[t 0 ,ti] C [0,T] and for all test functions ip £ C' 0,1 ([0, T] x Cl) 


0 = / / p(ipt + u T C7ip — c<p) + fpdxdt 

J tQ J 


pip(u T is)du)dt 


'to 


+ 


p(t 0 )g>(t 0 )dx 


p(ti)ip{h)dx 


( 10 ) 


holds. 
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Renormalization property: For any continuous and piecewise C 1 function /3 
the renormalization property holds, i.e. p satisfies for any [to>ii] C [0,T] and 
for all test functions <p G C' o, 1 ([0,T] x fi) 


0 = 


+ 


f f P{p)(Ft +u T V(p)dxdt - f f j3'(p) (pc- f)<pdxdt 
J to J ci J to J ci 


'to JCl 

f J 

'to J Cl 


ipdiv{u)(/3'(p)p — (3{p))dxdt — 


'to J r 


Pip/ p)(u T v)ipdu)dt 


/ P(p(t 0 ))<p(to)dx- / P(p(ti))ip(ti)dx. 
In J n 


( 11 ) 


Proof. See [3]. 


□ 


Remark 3. A component-by-component consideration shows the validity of Theorem [5] 
for vector-valued equations with n > 1. The renormalization property then reads: 

There exists a trace qp G L°°(Tt, |d/z u |) such that for any continuous and piecewise 
C 1 function /3: —> K." with D/3 diagonal, for any [to,ti] C [0,T] and for any test 

function tp G C , 0 , 1 ([0, T\ x 0,M") it holds 


0 = 


[ f /3(p) T (<p t +\3tpx)dxdt - f tp T (CD/3(p)p - D/3(p)f)dxdt 

J to J Cl J to 

f f ip T \J x (D/3(p)p — /3{p))dxdt — f f /3(7p) T (uU)<pdo;dt 

J to J Cl Jto JT 


/ /3(p(t 0 )) (fi(to)dx- / /3(p(ti)) T p(ti)da;. 
In Jn 


This theorem is no existence result, it only classifies solutions. But due to the 
continuity of p with values in L p (fl) and the existence of the trace we can define a 
solution of the initial-boundary-value problem as following: 


Definition 4. The function 


p G L°°((0, T) x fi) n 

is called a solution of the initial-boundary-value problem m if and only if 

- p is a weak solution of the transport equation, 

- the initial conditions are fulfilled, i.e. p( 0 ) = po and 

- the trace qp satisfies the boundary conditions, i.e. H{ r yp) = 7 p. 

As a consequence of the renormalization property we can prove the uniqueness of 
the above defined solution. 

Theorem 5 (Uniqueness). Under the assumptions from Sectionthere is at most 
one solution of the initial-boundary-value problem ©• 

Proof. Let p 1 and p 2 be two solutions of © with its traces qpi and 7 P 2 and define 
p = p 1 — p 2 . Then, p is a weak solution of the homogeneous transport equation. 
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Because of the uniqueness of the trace, 7 p = 7 P 1 — 7^2 is the trace of p. Thus, we see 
that p is a solution of the following homogeneous initial-boundary-value problem 

Pt + (U p) x + Cp = 0 in (0, T) x 

p( 0 , •) = 0 in fl 

(vU)-p = (v\J)-g(p\r T ) on re¬ 


using the renormalization property for (3(s) = |s| and ip = W1 with 1 = (1 
we conclude for all t. £ [0,T] 

0 = - / l T W/3(p(t))dx -II l T WC/3(p)da;dt 


Jn Jo Jn 

- I I l T W(uU)/3( 7/9 )dwdt 

Jo Jr 

= -||Wp(t)|| iW - f I l T WC|p|dxdf 
Jo Jn 

\ + 


I I |7P| T W( 
Jo Jr 


vuy — |^( 7 p)| W (vu) dcodt. 


With inequality ([ 6 ]) we get 




l|Wp(t)|| il(n) » < - 


0 J n 


l T WC|p|da;dt 


and thus 

l|Wp(i)|Ul (n) » < f \\c-\\ L oo {n) n\\W P \\ Llm ndt 

Jo 

< j (||(«X +c) _ ||i=o(p)n + ||(u x ) + || L oo(p)n^||Wp|| L l(n)ndt. 

By Gronwall’s inequality we conclude ||W / o(t)|| i i(Q)n = 0 for all t and thus p = 0 
almost everywhere since W is positive definite. □ 


Remark 6 . We can slightly weaken the assumptions on the source term /. In fact, 
for the existence of the trace, the renormalization property and the uniqueness it is 
sufficient to require / £ L 1 ((0,T) x Q) n (see [3]). 


Now, we are able to prove an upper and lower bound of the defined solution. 
Lemma 7 (Upper bound). Let the assumptions from Section^ be valid. Then, a 
solution p of the initial-boundary-value problem HD and its trace 7 p are componentwise 
bounded, i.e. it holds for all t £ [0,T] 


0 < p(t, •) < (p max H"~ F+{t)) ex P (^J a(s)ds^j 


almost everywhere in and 


a(s)ds 


( 12 ) 


0 < 7 p(t,u) < (p max + F+(t)) exp 


(13) 
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for \dp u \ -almost all ( t,u ) £ T^. Here, it is aft) 


(u x (t,-) + c(t,-)) 


and 

L°°( fi)» 


F+{t) 




/\\(h(s,-)) + h~\ 

OI(/»(*.0) + IIW 


Proof. Let ip £ C 0,1 ([0, T\ x fi,R”) be an arbitrary test function with p{T) = 0. In 
order to prove that rp with r{t) = exp J* a(s)ds^ solves also a transport equation, 
we would like to use ipr as a test function. Since r is only absolutely continuous, but 
not necessarily Lipschitz continuous, we need to approximate r by smoother functions. 
Therefore, let ak £ C' 1 ([0,T]) be a sequence converging to a in L 1 ((0,T)). We define 

afc(s)ds 

The mean value theorem states for all a<6el the existence of ( £ [a, b] with 

exp(6) — exp(a) = exp(£)(6 — a). 

Thus, we find for each t, £ [0, T] a constant f(t) with 

\C(t)\ <C k := max(||a fc || L i(( 0 ,T)), ||a|Ui((o,T))) 


r k (t) = exp (-/ 


such that it holds 


\r k (t)-r{t)\ = exp(C (t)) 


/ a k (s) — a(s)ds 
JO 


< exp(C k )\\a k - a||ii(( 0 ,T))- 

The sequence C k is bounded as ||afc|| i i(( 0i r)) is convergent, thus r k —> r in T°°((0, T)). 
Especially, it also holds a k r k —> ar in L 1 ((0,T)). Now, we will use < pr k as test 
functions in the weak formulation of the transport equation and take the limit for 
k —> oo. Because of ( r k )t = —oi k r k this yields 

0 = / / r k p T {<pt + \5<p x - (C + a fc Id) <p) + r k f T (pdxdt 

Jo Jn 

/ r k "/p T (vV)(pdu)dt + / pf)(p(0)dx 
Jo Jr Jn 

-A f f rp T (ip t + Up x — (C + aid) p) + rf T ipdxdt 
Jo Jn 

rT 


r"jp T (v\J)ipdudt + / pQip(0)dx 
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[ (rp-F + f (tpt + U<p, 
Jn 

T f (rf — (F + ) t ) T ip — F_ 
Jn 


(C + aid) ip) dxdt 
(U, + C + aid) fdxdt 


[ [ i r lP ~ F + ) T (vXd)fdudt + f p^p{d)dx. 

Jo Jr Jn 


Keeping in mind Remark^ this shows that p(t, x) = r(t)p(t, x) — F + (t ) with its trace 
r yp(t,uj) = r(t) 7 p(t) w) — F + (t) is the unique solution of the following initial-boundary- 
value problem: 


Pt + (Up), + Cp = f 

p(0,x) = p 0 (x) 

(vV)-p = (vU)-Hti>\r T ) 


in (0, T) x fl 
in SI 

dp“ on Tt 


with reaction term c = c+al, boundary operator T-L(p) = rp m — F+ + Q(F+) + Q{p) and 
source term /» = rf t - ||(r/j) + ||z,°°(n) - ((«»), + Cj) (-F+)*. For the boundary operator 
fi this conclusion needs further explanations. Since p is a solution of © we find 


7p|r T = np\r T - F + 

= rK(7p|r T ) - F+ 

= rpin + G(np\r T ) - F + 

= rpi n -F + + G{F + ) + G<dp\r T ) 
= JJ{lP\r T )- 


The advantage of introducing p is that proving the upper bound ED reduces to proving 
p < Pmax, where p max is known from the assumptions in Section [2] To this end, we 
define /3 : R" —> R” as 

/^(s) — (s Pmax)~^ 

and we will show /3(p) = 0 almost everywhere. 

For the test function ip = W1 and all to £ [0, T], the renormalization property 
yields 


/ " / 1 T W (CDp(p)p - DP(p)f) dxdt - / ° / l T WU x (D/3(p)p- p(p))dxdt 
Jo Jn Jo Jn 


/o Jn Jo Jn 

[ ( /3(7p) T W(Vu)dwdt — [ l T W/?(p(t 0 ))da; = 0 

/o Jr Jn 


(14) 


because of (3(po) = 0. The boundary term of this equation is non-positive. This can 
be shown using the componentwise monotonicity of /?, the assumption © and the 
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inequalities © and 
ft o r 


0{H( / yp)) W (yu) duidt 


/ o Jr 
to 


/3 ( rp in -F++ g(F+) + g( 7 p)) T W(u U )-dwdt 


< 


o Jr 

to 

o Jr 
rto 

lo Jr 

rto r 


/3 ( rp in - F + + G(F + ) + Q (p ma x) + G{lP ~ Pmax)) W( VU ) dwdf 


/3(Pmax + 0(7P-Pmax)) W (vu) dudt 


< 


< 


/o Jr 

, 

I o Jr 
rto 


W (yu) dcudt 


{^Q(lP PmaxJ 

G ({lP - Pmax) + ) W(vu)~diodt 

J /((7,--p m ax) + ) W(i/M) + dwdt 


/ o Jr 

= f j 0( r yp) T W(vu) + dujdt. 

Jo Jr 

Because of 0(s) > 0, D0(s)s > 0, c + u x = c + u. 

(■ D((3)f)i = D{f3)u (rfi - ||(r/ i ) + || i =o ( o) - {{ui) x + Ci)(F + )i) < 0 we conclude from 
equation © that it holds 

[ l T W0(p(t o ))dx = - / ° / 1 T W (C + U x ) D,3(p)pdxdt 


c + u x ) || L =c(n)nl > 0 and 


Jo J n 

rto 


rto 


l T WU x 0(p)dxdt 


< 


/ / VWD0(p)fdxdt + 

I o J n Jo Jn 

f f /3(7p) T W(i/u) + do;dt + f f !3(fi(^p)) T ^W(vu)~ dujdt 
/o Jr Jo Jr 

rto r 

Ti 


/o 


Wx) + ||i~(fi)” / 1 W/3(p)da;dt. 


/o 


(15) 


With Gronwall’s inequality this leads to f Q l T W0(p(t o ))dx = 0 and thus 0(p(to)) = 0 
and p(t 0 ) < Pmax almost everywhere. Using the renormalization property and the 
uniqueness of the trace of 0(p) = 0 , we also conclude 

0 = 7/3 (p) = /3 ( 7 p) 

and thus 7 p(t,w) < p m ax holds for \dp u |-almost all Addition of F and multipli¬ 

cation by exp ^ J* a(s)ds) shows the desired upper bound for p and yp. 

In order to prove the lower bound, we use the same procedure for p = rp and 
0 (s) = s~. In this case, it is c = c + a 1 , f t = rfi and 'H(p) = rp ln + G(p)- The 
only step we have to take into account again is the transition from equation (1141) to 
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inequality (ED- With the monotonicity of /?, the positivity of p m and the inequalities 
© and ©, it follows 

C *o . _ _ 

/3(H(7p)) T W (uw) _ dwdf 
/ /? (t Pin + GinP)) T w {vu)~dujdt 


J (jrp in + Ginp)^ ^ W(um) dwdf 


< 


< 


< 



Therefore, we can conclude as before p(to) > 0 almost everywhere in Tl and 7 p(t, oj) > 0 
for |dp. u |-almost all (t,oj) £ Tt■ This completes the proof. □ 

Beside the upper bound we just proved a lower bound. This lower bound can be 
sharpened if we ask for an additional assumption on Q and use the fact p > 0 almost 
everywhere. 

Lemma 8 (Lower bound). In addition to the requirements from Section © we as¬ 
sume ( c+u x ) £ L 1 ((0, T), L°°(Tl) n ). Furthermore, let p m i n £ K .™ 0 be a vector such 
that it holds 


(16) 


Pin exp C(s)ds^ +G(p min) ^ Pn 
component-by component for dp~ -almost all ( t,uj ) £ Tt and 

Po(x) > Pmin 

for almost all x £ Tl with £(s) = ||(« x (s, •) + c(s, •)) + lli oo (n)" ■ Then, any solution of 
the initial-boundary-value problem © is bounded from below. More precisely, it is 

rt 


almost everywhere in Tl and 


Pit, •) > Pmin ex p(-/ CO)ds^>0 
\nd 

Tp(t,u) > Pmin exp (^~ J q CO)ds^ > 0 


for all t £ [ 0 ,T] and |dp u | -almost all (f,w) £ Tt- 


Proof. Keeping in mind the non-negativity of p, this lemma can be proven in exactly 
the same manner as Lemma [7] □ 
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4. Continuous Dependence on the data 

The previous statements, especially the renormalization property of Theorem [5] 
provide us with the necessary tools to prove a central result of this article: a kind 
of sequential continuity of the solution operator. We formulate this statement as 
a theorem with two parts. In the first part, the weak-* convergence of a sequence 
of solutions to a solution of the limit problem is obtained for very weak assumptions. 
This result will play a crucial role in the proof of existence of solutions of the transport 
equation, whereas the second part, proving the convergence in C([0, T], L p (Cl) n ) under 
more restrictive assumptions, can for example be used for the proof of the existence 
of solutions of the low Mach number equations on a network (see » 

For the proof of the second part, we need the following result from the theory 
of Banach spaces concerning the uniform convergence of sequences of Banach valued 
functions. A proof can be found in El- 

Theorem 9. Let V be a uniformly convex and uniformly smooth Banach space. Let 
( f n ) n C C([0, T], Id) be a sequence and f £ C([0, T], Id) a function. Let the pair 
(( fn)n , f) fulfil the following two properties: 

(1) ||/n(i)||u converges uniformly to ||/(t)||y and 

(2) (tp(t), f n (t)) converges uniformly to (<p(t),f(t)} for all ip £ C([0, T], V'). 

Then, f n converges to f in the norm of C([0, T], V), i.e. 

lim sup \\f n (t) - f(t)\\v = 0. 

n-s-°o tg [o T ] 

Proof. See [T7]. □ 

Theorem 10 (Stability). 1. For all k £ N let Uk, Ck, fk, Po,k, Pin,k and Qk be defined 
as in Section 0 Assume there exists for each k a solution 

p k £ L°°((0,T) x fl) n 
of the initial-boundary-value problem 
( Pk)t T (U kPk\x T ^AkPk — fk 
Pk (0) = P0,k 

(uU k)~ Pk = (i / XJk)~'Hk(pk\r T ) 

with trace £ L°°(T t , |d/z„ fc |). 

Moreover, we assume: 

uniform boundedness: The sequence (p m ax, fc )fc C K> 0 from assumption (0 is 
bounded, i.e. there is p max £ R> 0 with p m ax,fc < Pmax- 
convergence of reaction terms: The sequence [ck)k strongly converges to c £ 
L 1 ((0,T) x f2) n in L 1 ((0,T) x f2) n . 

convergence of the velocities: (uk)k strongly converges in L 1 ((0,T) x f V) n to 
u £ L 1 ((0, T), Bd 1 ’ 1 (f2) n ) with (u x + c)~ £ L 1 ((0, T),L°°(Ll) n ) and (u x )+ £ 
L 1 ((0, T), L°°(fl) n ). 

convergence of velocities at the boundary: (vuk)k strongly converges to vu 
in L 1 {Y T ) n . 


in (0, T) x Ll 
in Q 
on Ft 
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boundedness of reaction terms and velocity derivatives: The sequence 
((cfc + ( Uk)x)~)k is bounded in L 1 ((0, T), L°°(Q)). 
weak convergence and boundedness of source terms: The sequence (fk)k 
weakly converges to f £ L 1 ((0, T),in L 1 (( 0,T) x Q) n and it is 
bounded in L 1 ((0,T), 

weak-* convergence of the initial conditions: The sequence ( Po,k)k weak¬ 
ly--*; converges to po £ L°°( fl) n in 

weak-* convergence of the inflow densities: The sequence (p ul .k)k weakly--* 
converges to Pin £ L°°(r T ,dp u ) in L°°(T T ,dp u ). 
weak-* convergence/continuity of the boundary operators: There exists 
an operator Q fulfilling the assumptions of Section^ such that for each weak-* 
convergent sequence qk £ -L°°(IY) n with limit q it holds 

C v\3k)~Gk{qk){vuk)~ (uU )~G(q) 


in L 1 (T T ) n . 

Then, there exists a solution p £ C([0, T\, L p (£l) n ) with trace 'yp £ |dp u |) of 

the transport equation 

p t + (U )p x + C p = f in (0, T) x 

p(0) = po in fl (17) 

(uU)-p = (vU)-H(p\r T ) onT T . 

Furthermore, it holds 


Pk P 

in L°°((0, T) x fl) n 

IPk A 7 P 

in L°°(T t , |dp„|) 

Pk(t) ->• p(t) 

in L p (n) n 


for all t £ [0, T] and p £ [1, oo). 

2. We assume additionally: 

convergence of the velocity derivative: The sequence (( Uk) x )k strongly con¬ 
verges to u x in the L 1 (( 0,T) x fi) ra -norm. 

convergence of reaction terms and velocity derivatives: (dk)k with ak = 
||(2 Ck + {uk)x)~ lli=o(n) n strongly converges to a = ||(2c + u x )~ Ili=o(o)n in the 
L 1 (0, T)-norm. 

convergence of source terms: The sequence {fk)k strongly converges to f in 
L 1 ((0,T) x fl) n . 

convergence of initial conditions: The sequence (po.k)k strongly converges to 
po in the L 1 (£l) n -norm. 
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weak lower semi-continuity of the boundary operators: For each weak-* 
convergent sequence pk —*• p in L 00 (^ 7 ’) r^ it holds 

f j /3(p) T W(uu) + — /3(%(p)) T W(z*u) _ dwdf 

j0 ^ rji (18) 


< liminf 


P(p k ) T W(vu k ) + - P(H k (p k )) T W(vu k )-dLodt 


with f3: R n —> R" defined by (3j(s) = sj, 'H k (p) = r k p- ln ^ k + G k (p) and 'H(p) = 
rpin + G{p)- Here, r k and r are given by r k (t) = exp 1 J Q t dfc(s)ds^ and 
r(t ) = exp 1 f* d(s)ds^ and W is the weight matrix introduced in Section [H 


Then, the convergence is even stronger: 
Pk t P 

and 


m C([0,T],L p (Cl) n ) 


TPk -> 7 P 

for all p £ [ 1 , oo). 


in L p (T t , \dp u \) 


Proof. As before, we extend ^p k , Pi n ,k, GkilPk) and Tlfi)p k ) by zero to functions in 

L°°(r r ) n . 

1. First, we want to prove the existence of a solution of the limit problem and the 
weak-* convergence. Because of Lemma [7] and the boundedness of p maXfc < /? max , 
II(c fc + (Mfc)x) _ ||L 1 ((o,T),L“(n)") < C x and ||/fc||n((o,T),i,«(n)») < C 2 the sequences p k 
and 7 p k are bounded almost everywhere by 

max (||pfc||i«>( ( 0 ,T )X fi ) ") ||7Pfe||L“(r T )-) 

< max ((p max ,fc + F k (t),F k (T) + F' fc (t))exp (||(c fc + {u k ) x )~ ||ii((o,T),L“(n)")) 

< (max(p max , C 2 ) + C 2 ) exp(C'i) (19) 

— Pm ax • 


Here, it is 


F k (t) = J exp j || (c k + (u k ) x ) ||dr 


/ ll((/fc)i(s,-)) + lli-(n) 


\ll((/fc)n(s, -)) + II ■) 


ds 


Thus, we can extract weak-* convergent subsequences - also denoted by p k and 7 p k - 
with 


Pk ^ p 

and 

7 Pk q 


in L°°((0,T) x Tt) n 


in L°°(T T ) n ■ 


As a second step, we will show that p is a solution of the transport equation and that 
q is its trace and fulfils the boundary conditions. 



TRANSPORT EQUATIONS ON A NETWORK 


16 


Because of the strong convergence of Uk and Ck in L 1 ((0,T) x f2) n , the sequences 
U kPk and CkPk converge weakly in L 1 ((0, T) x fi) n , i.e. 

U k pk U p 

and 


C kPk —*■ Cp. 


Similarly, (uUfc)ypfc converges weakly to (z/U )p in L 1 (Tt) u ■ Thus, for a test function 
<p £ C'°’^([0,T] x n) with ip(T) = 0 it holds 


0 = / pi {<pt + u kPx - C fc <p) + fk pdxdt 

Jo Jn 

+ / Po, k P>{ 0)d.T- / / 'ypl(vU k )ipdwdt 

Jn Jo Jr 

->• [ [ p T (<Pt+ U<p x - C<p) + f T pdxdt 

Jo Jn 

+ / pQ</j(0)da;— / / <7 T (i'U)<pdwdt. 

./n ./o ir 


In words, this means that p is a solution of the transport equation with its unique 
trace yp = q. Furthermore, because of Theorem [2] it is p £ C([0, T\, L p (tt) n ) and p 
fulfils the initial condition p(0) = po . For the boundary term we have to consider the 
trace q in more detail. The assumptions on Qk and pi n ,/c yield the weak convergence 
of (vUk)~Pin,k and (uU fc )~£fc(yp fe ). Thus, 



/ q (uU) _ <pdcudt t— / / yp fc (z/Ufc) _ <pdwdt 

Jr Jo Jr 

r T r 

'H(jpk) T {vUk)~vdujdt 


io Jr 
rT r 


/o Jr 

rT 


—> 


io Jr 

rT 


(Pin,k + QkilPk)) (uU k ) pdujdt 

j (Pin + G(,q)) T (uU)“(pdwdt 


'H(q) T (i' U) ipdudt 


/ o Jr 


is true and the trace q = yp fulfils the boundary condition 


OU) y q = (uU) Tli'yq). 


Hence, p is the unique solution of the limit problem (1171) . Since the weak-* limits p 
and q are unique, in fact the whole sequences and not only subsequences converge. 

As last step of the first statement, we have to show the weak convergence of Pk(t) 
for each t £ [0,T]. However, this is rather simple, as the sequence Pk(t) is bounded in 
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L p (f2) ra because of the continuity of p k with values in L p (JV) n and Lemma [TJ Thus, 
there exists a weak convergent subsequence in L p (Q) n with 

Pk(t) q- 

For this subsequence and for all time-independent test functions ip € C 0,1 (f2,R n ) it 
holds with Theorem [2] 

/ q T tpdx •<— / pk(t) T <pdx 

Jn J n 

n Pk (U kWx - C k p) + flipdxdt 
! 

+ / / / IPk (uUfc)vjdwdt 

Jn Jo Jr 

[ f P T (Utp x - C cp) + f T pdxdt 

Jo Jn 

+ / Pq pdx — / / 7 p T (uU)(/?do;dt 

Jn Jo Jr 

= / p(t) T <pdx. 

Jn 


Thus, we conclude 


q = p{t) 


due to the density of C 0 , 1 (f2, R n ) in L p (f 1) n . Again, by the uniqueness of the solution 
p we see that the whole sequence converges weakly. 

2. Now, we consider the stronger assumptions in order to prove both, the uniform 
convergence of p k with values in L p (fl)” and the strong convergence of the trace. As 
a first step, we prove that the convergence of Pk(t) is actually strong in L p (VL) n for 
each t G [0, T], To this end, we use the Radon-Riesz property (see e.g. [T3j) of the 
L p -spaces with p e (l,oo), i.e. the fact 


fk f in L p (n) n 

||/fe||Lp(n)" —► ||/||iJ>(n)" 


fk^f m L p (n) n . 


( 20 ) 


We define the auxiliary variable 

Pk(t,x) = Pk{t,x)r k {t) 

with r k {t) = exp (-5 f* d fe (s)ds^ and d k (s) = ( 2 c k + (u k ) x ) 


tion, it holds a k —>• a = 


(2c T iix) 


l°°( n)« 


L°°( Q) 

in T 1 ((0,T)). Thus, we conclude 


. By assump- 


1 


*(s)d,‘ 


rk(t) ->■ r(t) = exp 
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in C([0,T]), as in the proof of Lemma [7] We also see, as in that proof, that p k solves 
the equation 

(Pfc)t + (U kPk) x + C k pk = h in (0, T) xO 

jOfc(0) = p 0 ,k in fl 

(uUfe) _ pfc = (uU / t)^'Hfc(pfc|r T ) on r T 

with c fe = c fc + ia fc l, fk = r k fk and H k (p) = r k {t)p in ,k + Gk{p)- For p k and this 
problem, all requirements of the first part of this theorem are fulfilled. Thus, we 
conclude 

Pk(t) pit) 

in L p (fl) n , where p is the solution of 

(p) t + (Vp) x + Cp = I in (0, T) x fl 

p( o) = Po in fl 

{v\5)~p k = {v\J)-U{p\ TT ) onlY 

with c = c + |dl, / = rf and Hip) = r{t)p- m + Gip). 

The renormalization property for /3 : M n —>• 1” defined by /3(s)j = sj and = Wl, 
the weak lower semi-continuity of the L 2 -Norm, the strong convergence of po,fc and the 
weak convergence of /JW p k yield 


W/3(p 0 )dx — f f /3(7p) T W(i/M)do;dt 
Jo Jr 

(2c + u x ) T WPip)dxdt - 2l T Wpdxdt 
= I l r WPipit))dx 


/o Jn 


= l|Wp(i)|| L 2 (n)« 

< liminf ||Wydfc (*)||| 2 (rs) r, 


( 21 ) 


liminf l 1 W/3ip 0 , k )dx - J J Pi'yp k ) T W{u k v)dudt 


o Jn 


((2c fc + u k ) x ) T Wfiip k ) - 2 fl 


Wp k dxdt ] 


= / l J W/3( Po )dx 

Jn 

+ lim inf 

k \J o Jr 


/3(%(7Pfc)) T W(w fe u) - /3(7Pfe) T W iu k v) + dudt 


(2 Cfc + (u k ) x ) T W/3ip k )dxdt 


o Jn 




2 f 1 Wpdxdt 


o Jn 
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and thus 


limsupf [ [ /3(7p fe ) T W('Ufcu) + -/3(Hfe(7p fe )) T W(u fc u) dwdf 

k yJ o Jr 


+ 11 (2cfc + (uk)xY W f3(pk)dxdt 
/ o J n 


1 


( 22 ) 


< 


o Jr 


o Jr 


/ 9(7P) t W(uu) + - /3(W(7p)) T W(uu)"dwd< 

+ f f (2c +u x ) T W/3(p)dxdt. 

Jo Jn 

Using the assumption m for the weak-* convergent sequence X[o,t] 7 Pfc, we conclude 

/ 3{ r )p) T 'W(yu) + — /3('H( r yp)) T W(vu)~dujdt 

t (23) 

< liminf [ [ f3(pip k ) T W(vuk) + - / 3 (H fc ( 7 p fe )) T W(uu fe ) - dwdf 

k J o Jr 

for all t G [0, T). Because of the strong convergence of Ck and ( Uk) x in i 1 ((0, T) x fi) ra , 
we have for j = 1 ,..., n also the convergence 

24 + i u i)x 2 ~+ I 2 ? + (V)x | 2 

in L 2 ([0,T] x fl) and hence also the weak convergence 

24 + K)z 2 Pk | 2? + <y )*| 2 

in L 2 ((0,T) x fi). Here, the superscripts denote the components of the vectors. Since 
it holds 2 c k + ( u k )x = 2 c fc + (u fc )x + ||( 2 c fe + (u k ) x )~ ||L°°(n)” 1 > 0 , we can again use 
the weak lower semi-continuity of the norm to find 

pt p n 

/ / (2c + ^fW/3(p)dxdt = ^||(2? + K) x )5W^|| 2 2((M 
J o J r 2 j 

n 

< limhrf^ ||(24 + K)x)^W ii 4 ||| 2((0it)xn) 

3 =1 
r t 


= liminf / / (2c fc + (u fc ) x ) W/j(p fe )da;dt. 


o Jn 


This leads with the inequalities (E2> and (E31) to 


lim sup 

k 


. , , /3(7 Pfe) T W(u fc u) + -/3 (W fe (7/3 fc )) T W(u fc u) dwdf 

\do dr 

(2c fe + (u fc ) a: ) T W/3(p fe )da;d< ] 
o Jn / 
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< 


/3(7 p) W(uu) + — (3('H( , yp)) W(uit) dwelt 


o Jr 


+ / (2 c + u x ) W/3(p)da;df 

J o J n 


< lim inf 

k \J o Jr 


PilPk) W(u fc u)+ - Pi'HkhPk)) W(M fe u) dw dl 


+ lim inf / / (2c fc + (-Ufe) x ) T W^(/9 fe )da:dt 


< lim inf ^ 


/o ./ft 
r t 


Php k ) T W(u k v) + - mkilPk)) 1 VJ{u k v)~dujdt 


- 


o ./r 


+ J J (2c fc + (Mfe)x) T W/3(p fc )d.Tdf j . 

Hence, the limit inferior and the limit superior coincide and consequently they are 
equal to the limit. Altogether, this yields 


k ^o ll w ^)lllw 


= lim - 
k—yoo \ 


/3(7Pfc) Tw («feu) + -^(Hfc(7p fc )) T W(u fc i/) dwdt 


o 


l T W/3(p 0 , fc )dx - 


(2cfe + (u k ) x y Wp(p k ) - 2/^ W/ifedxdt 


o Jn 


/3(7 P) j W(uu) + - Pi'Hi'yp)) 1 W(vu)~dujdt 

o Jr 

+ [ l T W/3(p 0 )dx- [ [ (2c + u x ) T W/3(p) -2f T Wpdxdt 
JO JO JO 

= l|Wp(t)||| 2(n) „. 

Now, the strong convergence of p k [t) for all t £ [0,T] follows by the Radon-Riesz 
property (1201) of L 2 (H) r^ and the equivalence of the weighted Euclidean norm ||W • ((2 
and the Euclidean norm || ■ || 2 - Thus, by the convergence of r k 

r k (t) = exp Q J d k (s)ds' S j -A f(t) = exp Q J a(s)di 

in C([0,T]) (see proof of Lemma 0, the convergence 

p k (t) = p k (t)f k (t) -A- p{t)r(t) = p(t) 


in L 2 (Q) n is also true. Because of the uniform boundedness of p k {t) in L°°(Q ) n , we 
see that p k (t) converges in all L p (B) n -spaces for p £ [1, 00 ). 
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Up to now, we have shown the pointwise in time convergence of pk■ In the next 
step, we want to show that the convergence of pk is uniform in [ 0 ,T] with values in 
L p {Ci) n . This can be done by applying Theorem [9] to pk, a Radon-Riesz like property 
for uniform convergence. Therefore, we first prove that f n p k (t, x) T <p(t, x)dx converges 
uniformly towards f Q p(t, x) T ip(t, x)dx for all fixed tp G C' 1 ([0,T] x Recall the 

upper bound p max G R>o for pk from the first part of the proof and let be £ > 0. We 
define £ = — n—n— -• Since Uk, Ck , fk and vu k are converging in L 1 , they are 

0Pmax|| ( P||c' 1 ([0,T] xfi,K n ) 

equi-integrable. Thus, there exists 6 > 0 such that S C [0, T\ with p(S') < 5 implies the 
four inequalities: f s ||ufc||z,i(n)ndt < £, f s ||c fc ||ii(n)™dt < £, f g ||/fc||.Li(fi)«df < £ and 
f s ||uufc||ii(r)'idt < £. Then, we compute for ti,t 2 G [0,T] with \t 2 — ti| < min(<5,e) 


p(t 2 ) T p k (t 2 )dx - / <p(ti) T p fc (ti)dx 
J n 

't-z p pi2 p 

/ (pf + V?zUfc - p T C k ) Pk + flpdxdt - / / p T (vU k )"/p k dujdt 

i i/fi i/•/r 


< IMIcrPmax / (l + llwfclUqn)" + ||cfc|Ui(n)n + H/fclUqn)") df 

Jti 

pt2 

+ IMI ooPmax / ||uwfc||Li(r)«df 
Jti 
pt 2 

< PmaxIMIc 1 / (l + ll w fc||.L 1 (fi) n + || c k || L 1 (Q) 11 + ll/fc II L 1 (fi)" + \Wu k ||i 1 (r) n ) dt 

Jt i 


< £ 


which proves the weak equicontinuity of p k ■ Together with the weak convergence of 
Pk{t) in L p {fl), we conclude the uniform in time convergence of J n pk(t) T p(t)dx as 
each equicontinuous pointwise convergent sequence is uniformly convergent. Due to 
the density of C 1 ([0,T1 x SI)™ in U([0, T],L 2 (Q) n ), this uniform convergence holds in 
fact for all ip G C([0, T], L 2 (Q) n ). 

To apply Theorem [9j we additionally need the uniform convergence of the norm 
l|Pfc(OIU 2 (fi)"- Therefore, let again be e > 0 , define £ = 4max (p max)1) p max and choose 
a value S > 0 such that for all S C [0,T] with p(S) < S the inequalities f s \\2c k + 
(wfe)x||z,i(f2) ndi < fo Is WfkW^iQ^dt < £ and / s ||i/u fe || L i (r )ndt < £ hold. Then, the 
renormalization property leads to 

l|pfc(^ 2 )||| 2 (n)n — ||pfc(U)||i 2 (n)" 

= f f P{lPk) T {vu k )dudt + f f (2c k + {u k )x) T PiPk) ~2.fipkdxdt 
j t\ J r Jt\ J 

pt 2 

< max(p max , l)p m ax / (II VUk II Li(r) rl T||2Cfc + (lifc)a; T 2\\ fk\\L 1 (Q) n ) d t 

Jt 1 

< £ 


forti,t 2 G [0, T]with \t 2 — 1\ | <6. Again, the shown equicontinuity yields together with 
the pointwise convergence of ||jOfc(t)||L 2 (f 2 ) n the uniform convergence. As a consequence, 
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we can apply Theorem [9] to deduce the convergence of pk in C([0, T], L 2 (Q) n ). As 
before, the convergence holds in fact in all the spaces C([0, T], L p (fl) n ) for p £ [1, oo) 
because of the uniform boundedness. 

To finish the proof, it remains to show the strong convergence of the trace 7 p k - 
Again, we consider the renormalization property for /3 as before with an arbitrary 
<p £ C°’ 1 ([0,T] x Since fi(p k ) converges strongly to /3 ip) in L p ((0,T) x fl) n 

and (Kp k (T)) to /3(p(T)) in we conclude 


/3{lPk) T (vU k )<pdujdt 


/ 0 Jr 

r T 


f f /3{pk) T {<Pt + Vkfx - (2C fc + (U fc ) x )<p) + 2f£p k dxdt 
Jo Jn 

+ [ /3(p 0ik ) T <fi(0)dx - [ f3(p k (T)) T ip(T)dx 

Jn Jn 

f [ P{p) T (<p t + \]p x - (2C + U,V) + 2f T pdxdt 
JO Jn 

+ [ /3(po) T ^(0)dx - [ /3(p(T)) T ip(T)dx 

J n J n 


(24) 


/3 ( 7 / 3 ) T (uU) ipdwdt. 


>0 Jr 


To use the Radon-Riesz-property again, we would like to choose sgn(uzt) as a test 
function in order to receive a norm of yp. But since this function is not smooth 
enough, we have to approximate it to ensure the C 0,1 -regularity. Therefore, let be 
tfij £ C 0,1 ([0, T] x such that 


(fij -A sgn(uzi) 


holds almost everywhere in Ty with | ip(t, ir)| < 1. The dominated convergence theorem 
immediately yields 

(vXJ)tpj —> |uzt| 

in L 1 (T T ) n . Thus, we find for e > 0 an index J £ N with 

£ 

||(uU)<pj - \uu\\\ L i {rT) n < ——. 

max 

With the convergence in (liHl) we find for pj an index K\ such that 


[ [ (/3(7Pfc) T (^Ufc) - /3(7p) T (uU)) ipjdudt 

Jo Jr 



for all k > K\. Furthermore, because of the convergence of vu k in .^(Tt)" there is 
also K 2 with 


v{u k 


w)|Ui(r T )n 


£ 

4 Ip 

r max 
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for k > A'2. Finally, for all k > max(A'i, A'2) it holds 
ll7Pfc|lL 2 ( rr ,| d ft u |) “ ll7PllL2(r T ,|d/i u |) 

/ (3( r ypk) T {vU) (sgn (vu) - <pj) dwdf 


< 


0 Jr 

n 

/o Jr 

n 

/o Jr 

r T 


P(lPk) T (v (U - U fe ))<pjdwdt 
(P{lPk) T {u\Jk) ~ P{ip) T (uIJ)) (fjdudt 
/3(7p) T (j/U) (ipj - sgn(uzi)) dwdt 


0 Jr 

< 2 PmaxlllH - (^U)v?j|| L i ( r T )n + p 2 ma Ju(u - u k )\\ L i { r T r 

< £ 


and we can use the Radon-Riesz-property and the boundedness of 'ypk to conclude the 
strong convergence in L p (Tt, |d/i„|) for all p G [1, 00 ), which finishes the proof. □ 

Remark 11. Boyer and Fabrie obtained in [3] an equivalent result for the scalar trans¬ 
port equation with non-coupled boundary conditions. However, their proof technique 
is completely different and uses a mollifying procedure. 


5. Existence 

Now, we are finally able to prove the existence of solutions of the transport equation 
©. This will be done in two steps. In a first step, we will consider the equation with 
non-coupled boundary conditions. In this case, the existence of a solution is classical 
and it can be proven in several ways. One possibility is to use the method of char¬ 
acteristics for smooth data. Then, the first part of the previous theorem yields the 
existence for general data. In [3] , a proof of the existence using a parabolic approxima¬ 
tion is presented. In a second step, we will construct a solution for general boundary 
conditions by an iterative method, where we will use the first part of Theorem 1101 
Lemma 12 (Existence). For Q = 0 there exists a solution of the initial-boundary- 
value problem © under the assumptions from Section [3 

Proof. See e.g. [3] or m- n 

Now that we know the existence of solutions in the case of non-coupled boundary 
conditions, we will provide a small lemma concerning the monotonicity of the solution 
operator in dependence on the inflow values. This lemma will turn out to be useful for 
the construction of the solution in the general case with coupled boundary conditions. 
Lemma 13. ForQ = 0 let the assumptions from Section\^be true and let the functions 
Pin,l; Pin, 2 G A°°(Tt, d/i“) be two different inflow values with 

Pin,2 > Pin,l > 0 
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d p u -almost everywhere. Then, for the solutions pi of the transport equation © with 
inflow values (z/U) - /^ = Pi n ,i it holds 

P2(t) > Pl(t) 

almost everywhere in Cl. Moreover, for their traces jpi hold 

1P2 > IPi 

\dp u |-almost everywhere inTr- 

Proof. Let p± and P 2 be the solutions corresponding to the inflow values pi n! i and Pi n , 2 - 
Recall that these solutions p\ and P 2 exist due to Lemma [T2] Then, the difference 
p = P 2 — pi with trace 7 p = 7/32 — 7Pi solves the homogeneous transport equation 

Pt + (U p) x + Cp = 0 in (0, T) x Cl 

p( 0, •) = 0 in Cl 

(uU)“p = (uU)“ (pin,2 — Pin,l) OnTj-. 

Since the inflow value Pi n ,2 — Pin.i is non-negative dpf almost everywhere, we can apply 
Lemma [7] to obtain the lower bound p{t) > 0 almost everywhere in Cl and 7 p(t,u>) > 0 
for |d^„|-almost all (t,oj) £ Tt- □ 

Equipped with this lemma, we can prove the existence of solutions for general 
boundary operators TL using an iterative construction in the next theorem. 

Theorem 14 (Existence). Under the assumptions m on u, c and f and under the 
assumptions from Section [H on H, there exists a solution of the initial-boundary-value 
problem ©■ 

Proof. We will iteratively construct a sequence of functions converging weakly-* to a 
solution p of the transport equation. The idea is to use the outflow value of a solution 
to compute the inflow value for the problem of the next iteration step. Due to the 
positivity of the operator Q, the obtained sequence will be monotonically increasing. 

Therefore, we define the starting value of the sequence 7/30 = 0. From Lemma fl2l 
we know that the problem 

(Pk)t + (U p k )x + C p k = f in (0, T) xO 

Pk( 0) = Po in O 

p k = {v\3)~UflPk-i) on T t 

has a unique solution pk with trace 'ypk- Once more, we want to use the stability 
theorem ('Theorem flOl) . but this time with pi nt k = TLipipk-i) and Qk = 0. 

As a first step, we have to check the requirements for this theorem. Therefore, we 
want to analyse the boundedness of the boundary value. From assumption © we 
know that there is a constant p max G K” 0 with 

a(s)ds 




Pin exp 


F + G(F) + Q (/3 m ax) < Pmax 
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and 


Po - Pmax 


almost everywhere. We want to show by induction that an analogue bound for the 
boundary value = T~L{lPk-i) and Gk = 0 is valid for all k , i.e. 


a(s)ds 

The base case is clear since Q and p ln are positive and it holds ypo = 0. Using the 
positivity of Q, the induction hypothesis and the assumption fl£]) on p max we calculate 
for the inductive step 


^ / Pmax- 


7 Pk exp 


'H{lPk-i) exp J a(s)ds^ 

= Pin exp J a(s)ds' S j + G(p max ) + G(F) 

-G ^Pmax +-F 1 - 7Pfc-i exp J a(s)ds^ 
<Pi n e*p(^-J a(s)ds^ + G(p max ) + G(F) 


^ Pmax + F '. 


Thus, Lemma [7] yields the desired bound for 7 pk for all k. 

The next aim is to show the convergence of the inflow values pi n ,fc- As we have just 
seen, this sequence is bounded and it is monotonically increasing, i.e. it holds 


Pin,fc+1 7 Pin,k (t,w) 

for dp“-almost all ( t,uj ) £ Tt- We will also prove this by induction. Due to the 
non-negativity of / and Lemma 0 it is 7 pi(t,w) > 0 for |dp u |-almost all ( t,to ) and 
thus, because of the positivity of G and p; n 

Pin,2 = %(7Pl) > Pin = H(7 Po) = Pin.l 

holds dp“-almost everywhere. For the inductive step, we will use the previous lemma 
about the monotonicity. Assume, it is pi n ,fc > pi n ,fc— 1 - Then, Lemma [T3l immediately 
yields 7 pk > jpk-i- Thus, the positivity of G leads to 

Pin,fe+1 - Pin, A: = %( 7Pfe ) - 'H(”/p k -l) = G{lPk ~ IPk-l) > 0. 

Since (p; n ,fc)fe is bounded and monotone, we can apply the monotone convergence 
theorem to conclude that there exists a function q £ L°°(r 7 ’,dp“) with 

Pin ,k ^ Q 
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in L p (Tt, dp“). Since the sequence is bounded in L°°(Ft, dp“), it is also weakly-* 
convergent in L°° (Ft , dp~). Hence, Theorem Hill is applicable and it yields the exis¬ 
tence of a solution p with trace qp of 

p t + (Up)* + C p = f in (0, T) xfl 

p(0) = po in H 

(uU)“p = (uU) - g on T-r- 

Furthermore, it holds 

7 Pk 7P 

in L°°(r T , |dp.„|). The weak-* continuity of the operator Q implies 

ft(7Pfc) A %(7P) 

in L°°(r T ,dp-). Finally, we conclude 

(z/U) _ 7p (uU)“7p fc 

= (uU)-^( 7 p fc -i) 

- (vV)-H('yp) 

in L 1 (Ft'). Because of the uniqueness of the weak limit the boundary condition 

(uU)“7p = (uU)-%(7 p) 

is satisfied. □ 


Remark 15. Remark, that the required non-negativity of / is essential in the pre¬ 
vious proof. In order to prove the existence without this restriction, one needs an 
additional requirement on Q, e.g. the existence of a constant vector p m ; n G R" 0 with 
G{F- + p m in) < F_ + p min for dp~ almost all (f, w) G r T - Here, it is F_(t) = 


exp 


(- fo a(s)ds) 


ds. Then, we can split the problem in the two 


VII (fn(s,-)Y 


subproblems 


Pt + (Up 1 )* + Cp 1 = /+ 

in (0, T) x 

p x (0, ■) — po 

in Id 

(uU)-p 1 = (uU)'7£(p 1 |r T ) 

on Tt 

pt + (Up 2 )* + cp 2 = r 

in (0, T) x H 

p 2 ( o,-) = o 

in H 

(uU)-p 2 = (uU)-0(p 2 |r T ) 

on 


so that p = p 1 — p 2 is a solution of the original problem. Due to the existence of p m i n , 
both subproblems fulfill all requirements of section [2] 
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6. Application on networks 


One of the main applications of such coupled boundary value problems are transport 
problems on networks or graphs. In this case, the boundary operator describes the 
coupling of the density at the nodes of the network. As usual, we enforce the mass 
to be conserved at the nodes. Furthermore, the density should be distributed by a 


prescribed ratio. This is similar to the coupling conditions for traffic flow (see m or 


open channel flow (see [6] or p~5]l. 

In order to describe such flows on a network, the notation of graph theory is very 
well suited. We will briefly introduce the notions which are needed in this context. 
For more details, we refer to the PhD thesis m and to any text book about graph 
theory as e.g. [7|. 

Let G = (V, E, w, init, ter) be an oriented, weighted and connected graph with n 
edges and m nodes and let B £ { — 1,0, l} mxn denote its incidence matrix. Assume 
that the graph has k > 0 inner nodes Vj £ V with degree d(vj) > 1. We denote by 
B>i £ { —1,0, l} fcxra the submatrix of the incidence matrix containing only the rows 
of B corresponding to inner nodes. In the same way, B = i £ { — 1,0, i|( m - fc ) xn i s the 
submatrix corresponding to the outer nodes Vj £ V with degree d(vj) = 1. 

With this notation, we are able to write the mass conservation at the inner nodes 
as 



If there is more than one outflow edge at the same node, we need to specify the 
distribution of the density. The simplest model is the perfect mixing of the density, 
i.e. the density in the outflow edges is continuous at the node. However, also other 
linear models to determine the distribution on the different outflow edges are possible. 

In order to define the above-described boundary operator for the perfect mixing 
precisely, we introduce the matrix M £ L 1 ((0,T)) kxk by 



(25) 


A simple computation using the structure of the incidence matrix shows that this 
matrix M = ( niij ) has a diagonal form with 



(see US])- For the well-posedness of the transport equation we need additionally some 
sort of energy conservation. We assume 



(26) 


for all inner nodes v% £ V with d(vi) >1. In the formal low Mach number limit of the 
Euler or Navier-Stokes equation on a network studied in |16l this condition is fulfilled 
naturally. 
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Then, for p out G L°°((0,T)) m k with p out > p m i n > 0 almost everywhere, we define 
the boundary operator 

n u : J L 00 (r T ,d/ i +)^L 00 (r T ,d/ i -) 

P ^ Pin + 

with 


Pin = (pB^) Pout 


and 


<?«: I/°°(rT,dpu) —»• L^^Tjdp”) 

S.(p) = Kt) - M- 1 J (i/B >1 ) _ (z/U) + pdw. 

Here and in the following, MA 1 has to be understood in the sense of the Moore-Penrose 
pseudoinverse. Clearly, this boundary operator satisfies the assumption on the perfect 
mixing since for each inner node there is a common outflow value for the density. 
Furthermore, the mass is conserved at the inner node since it holds 


btdpjdw 


jf (uB >1 )- (uU)- Hu(p)du = j (vU)' 

= MM- 1 / (uB>i)" (uU) + pdw 
= ^(uB > 1 )-(uU) + pdw 


(27) 


for all p G T°°(r T , dp+) and almost all t G [0, T\. For the first equality we used the fact 
(uiB>i)-W(u 2 B; i )- = 0 for arbitrary diagonal matrices W and ui,u 2 G {—1,1}. 
And for the last equality we have to take into account the energy conservation (ITHl) 
since M- 1 is just the pseudoinverse: In the case m„; = 0 it is also 


J ^(uB>i) (vu) + ^j dw = 0 


by ([26]). This is a sum of positive terms and thus all summands have to be zero and 
the third equality in (1271) is obvious. 

Until now we have formally seen that the operator describes the above-mentioned 
coupling conditions. Before we check that H u fulfils the requirements of Section [2] for 
the well-posedness of the problem, we will show the well-definedness of Q u and we will 
introduce the pre-adjoint operator Q u , to which Q u is the adjoint operator. 

For all p G L°°(Tt, dp+) and for almost all t G [0, T] it is either rrijj(t ) = 0 or 


m jj(t ) 1 / f(uB>i) (z/U(t)) + p(t, w)') dw < max(pj(t, w)). 
Jr ^ '5 

In the case rrijj(t ) = 0, it is (Q u ) i (p)(t,w) = 0 for all adjacent edges e*. 
Altogether, this proves 


(28) 
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for almost all (f,cc) G Ty with vui(t,u>) < 0. Hence, it is Gu(p) G L 00 )^, d/z“) and 
the linear operator Q u is well-defined. Descriptively, this is not surprising since the 
coupling condition chooses a weighted mean value of the inflow densities for each t. 
Clearly, this weighted mean value is smaller than the maximum. 

As a next step, we define the operator Q u as 

Gu ■ L 1 (T T ,dp~) L 1 (T T ,dp+) 


p i-> M _1 y (uB>i) _ (uU) _ pdw. 


This operator is well-defined by the same argument. We observe that Gu is the adjoint 
operator of Gu- To check this, let be p € L°°(TTApt) and P G T 1 (Tt, d/x“). Then, 
we compute 

[ p T (vU)-G u (p)dcudt= [ [ p T {v\l)~ (uB^'docM" 1 f (uU) + pdu; 

Jy t J o Jr Jr 

= [ Qu{fi) T (^U) + pdudt. 

Jr T 

Now, we will check all requirements for the boundary operator assumed in Section [2] 
Lemma 16 (Continuity equation on a network). Let u € T 1 ((0, T), W 1,1 )!))") 
with u x G L 1 ((0, T), L°°(n) n ) and po G L°°(Vl) n he given and let be H u the above- 
defined operator. Then, the continuity equation on a network 

p t + (U p) x = 0 in (0, T) x Q. 

p{ 0, •) = po in fi 

(i/U)"p = (vUrH u (p\r T ) onV T 

has a unique solution which is bounded and depends continuously on po and p ln . In 
addition, for each po G L°°(H) n with po > p m j n > 0 almost everywhere, there is 
a vector p min G R" 0 satisfying inequality (fTB|) . which guarantees the existence of a 
strictly positive lower bound of the solution. 

Proof. We will verify each of the assumptions of Section [2] so that we can apply the 
theory of the previous sections. Therefore, let be p G L°°(rT,d^+). 

Weak-* continuity: Since Gu is the adjoint of another operator, namely Q v • it is 
weak-* continuous. 

L 1 -operator norm: As weight matrix for the L p (Tt, dp,±) norm we choose the 
identity matrix W = I. Then, using the definition of M, we estimate 


\\Gu{p)\\i, w - = [ \Gu(p)\ T (vu) dudt 

JFt 

p T (nV)+ (uB^) 


IYt 

r T 


p T (v U)+ (uB t 


>1/ 


doc 


doc 


M 


ldf 


(uB>i) (vu) docdf 


^ ||P||. 
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Thus, the .^-operator norm of Q u is less or equal one. 

Causality: Equation Q is fulfilled for almost all t £ [0,T] since it holds 

X{o,t]Gu{X{o,t]p) = M _1 J (uB>i)~ [v U) + X[o,t]pduj = X[o,t]G u (p) 

by definition of Q u . 

Positivity: Let be p > 0 almost everywhere. Then, Gu{p) is positive as sum of 
products of positive factors. 

Boundedness: For any po L°°(Gl) n define the scalar 

Pmax = max(||p out r|| L oo (rT) m-fc, ||p 0 ||i<»(fi)») 

with r = exp J 0 * ||(u x (s, -)) _ ll-L°°(fi)«ds^. Observe that the reaction term c and the 
source term / are not present in the continuity equation. Due to the definition of the 
matrix M, it is 

Gu(Pmax^-rt) — PmaxGu(ln) = Pmax (uB^) 1VE IVElfc = Pmax (uB^) lfc 
and thus, it follows also 

PinT "F ^u(Pinaxln) = PoutT “F Gu (Pmax^n) 

Ff Pmax (uB_^) lm-fc “F Pmax (uB^) lfc 

= Pmax(uB ) l m 

— Pmax In 


dp“-almost everywhere. Here, one has to be careful with the different dimensions of 
the 1-vectors. In the last step, we used that (z/B) _ has in each column exactly one 
non-vanishing entry. 

Now, we additionally assume po > p m in > 0 almost everywhere and denote 
f{t) = exp ||(w,(s,-)) + ||ds 

Then, we define 


Pmin = min{essinfp 0 ,essinf(rpout)} > Pmin > 0 
£7 r T 

and thus, it holds 


PinT + 5n(Pminlra) — (uB^) PouiT “F PminGu(ln) 

fF Pmin (uB_^) 1 m—k “F Pmin (uB^j) Ifc 

= Pmin(uB ) l m 
— Pmin In 


dp u -almost everywhere. This completes the proof. 

□ 


Remark 17. Using the stability theorem[lOl one can also proof that the solution of the 
equations on a network continuously depends on the velocity u. 
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